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We quantitatively investigate the quark mass dependence of current matrix elements and energies, calculated 
over a wide range of quark masses in the continuum limit of small- volume quenched lattice QCD. By a precise 
comparison of these observables as functions of the heavy quark mass with the predictions of HQET we are able 
to verify that their large quark mass behaviour is described by the effective theory. 



1. MOTIVATION 

The simplification of the QCD dynamics in the 
limit of large masses of c- and b-quarks gives rise 
to the Heavy Quark Effective Theory (HQET) 
as a standard phcnomcnological tool to describe 
decays of heavy-light hadrons and their transi- 
tions in terms of hadronic matrix elements. (See 
e.g. ^ for a review.) Starting from the HQET 
Lagrangian of a heavy quark, £hqet , that reads 
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this effective theory provides an expansion of the 
QCD amplitudes in the inverse heavy quark mass, 
1/m, and is renormalizable at any finite order in 
1 /m according to power counting. 

Making HQET an effective theory for QCD re- 
quires matching calculations to express the pa- 
rameters m, cjgpin, . . . in 1^ by those of QCD, 
and the agreement of different determinations of 
quantities such as Vch ^2 , which involve pertur- 
bative HQET, already reflect the success of the 
effective theory approach. Apart from these phe- 
nomenological tests of HQET, however, indepen- 
dent, non-perturbative ones are still desirable. 

Though in principle the lattice achieves this by 
varying m, clean comparisons of QCD and HQET 
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in the continuum limit demand m <^ 1/a prior to 
a — *■ 0. Yet j3] rather turned that restriction into 
an idea to connect them non-perturbatively: con- 
sider QCD in a small volume L^, where the b can 
be simulated as relativistic fermion and HQET 
becomes an expansion of QCD in the variable 

1/z = 1/{ML) , M: RGI heavy quark mass . 

Here, we confront the large-z behaviour of cor- 
relation functions, computed with Schrodinger 
Functional boundary conditions ^4,, to the static 
theory and estimate the size of 1/m-corrections. 
This is also of practical relevance in the strat- 
egy to solve renormalization problems in HQET 
non-perturbatively by matching to QCD in finite 
volume T . For a full report on our study (and a 
list of references to related work) consult jSj. 

2. LARGE-MASS ASYMPTOTICS 

Our observables derive from relativistic heavy- 
light SF correlation functions. With /a a correla- 
tor between a heavy-light pseudoscalar boundary 
source and the axial current Aq = '017075 V'b in 
the bulk, ky its vector channel analogue and /i a 
boundary-to-boundary correlation, we define 
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which from the multiplicative renormalizability of 
the SF boundary fields follow to be finite quan- 



1 



2 



titles provided that , Vq denote renormalized 
currents. Effective energies are constructed as 

rps(i, M) = -/;(r/2)//A(T/2) , ry: /a ^ fcy 

with spin-averaged sum and difference 

Fav = i[rps + 3rv] , Ar = rps-rv. (3) 

Being (ratios of) matrix elements between low- 
energy heavy- light and vacuum- like states 
© and Q should then be described by HQET, 
which we test by verifying their large-z asymp- 
totics to comply with the predictions of HQET. 

To obtain them, note that classically one ex- 
pects the current matrix elements to be power se- 
ries in 1/2; led by the static limit, where the heavy 
quark does not propagate in space. E.g., replac- 
ing in ^ ■0b by iph and dropping the 0(l/m) of 
'Chqet in Q]) to evaluate the static correlators, 
the effective theory and QCD are related by 



X{L) ^ /r*(i)/y7f^ = lim YpsiL, M) . (4) 

z — ^oo 

On the quantum level, the scale dependent renor- 
malization of the effective theory implies logarith- 
mic modifications, i.e. the axial current renormal- 
ization amounts X]x{L,^i) = Z^^'^^{ij)X\^s,re[L) to 
depend logarithmically on the chosen renormal- 
ization scale (/i) and scheme, but not so the as- 
sociated renormalization group invariant 
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[L)^ lim f[26og'(Ai)]-^XR(L,M)}, (5) 



calculable in lattice QCD: Xrgi = ^RGi-^bare IE]' 
The large-z behaviour of Q thus splits into RGIs 
of the effective theory and logarithmically mass 
dependent functions, C. As their arguments we 
choose r = M/Aj^, since it can be fixed on the 
lattice without perturbative uncertainties I7l8j : 

yps(£, M)^'^°^ Cpsir) Xrgi{L) + 0(l/z) , (6) 
R{L, M) ''^^ Cps/v{r) (1 + 0(l/z)) . (7) 
Similar predictions hold for the energies lO, 

Lrav(L, M) Cm^sir) z + 0(4) , (8) 
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where Cmass translates the pole into the RGI 
quark mass and A^rq" is a RGI matrix element 
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Figure 1. Continuum limit extrapolations span- 
ning the entire z-range. (Dotted lines extend the 
linear fits in {ajhY to omitted values of a/L.) 

of the chromomagnetic operator ipY^a-'BiJjh, whose 
anomalous dimension (AD) contributes to Cgpin- 
Taking the entering ADs to best perturbative 
accuracy [H], the conversion functions C in ©- 
are evaluated by solving RGEs |S] . As detailed 
there, their perturbative knowledge is adequate to 
allow for an investigation of the 1 / z"-corrections. 

3. DISCUSSION OF THE RESULTS 

Our quenched data refer to a volume of extent 
L = T « 0.2 fm, which admits to reach M > 2Mb 
while a — > extrapolations are still well control- 
lable. To account for the growing of the quark 
mass in lattice units at given a/L as z is increased, 
the coarsest resolutions that pass into the contin- 
uum extrapolations linear in {a/Vf' are chosen 
by imposing a cut on aM 5 8 . Fig. ^ illustrates 
a typical case: the slopes are quite small, and the 
error in the continuum limit gets the larger the 
more a/L are to be discarded for increasing z. 

The main results and their polynomial fits in 
\/z to quantify the deviations from the static 
limit are displayed in figs. [2-01 linear fits use only 
the heavier quark mass points. 

3.1. Current matrix elements 

Comparing the finite-mass matrix element of 
^0 with the HQET prediction Argi(^) + 0(l/z), 
we infer from fig. [21 that the perturbative Cps 
reduces the mass dependence of Yps significantly 
and renders Ypq{L , M) / Cps{f) to be compatible 
with approaching the static result for Arqi |5I6| 
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Figure 2. Fits ofYps/Cps include Xrgi- The AD 
of the static axial current (7 ) enters in Cps • 
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Figure 3. Fits of R/Cps/y, constrained to 1. 

as 1/z — > 0. Also the ratio R of matrix elements 
of Aq and Vq in fig.|3|is consistent with the leading 
term in the 1/z-expansion (fixed to 1 by the spin 
symmetry of the static theory, cf. (0), if Cps/v 
is evaluated including at least the current's two- 
loop ADs. In both cases, the coefficients of 1/z- 
corrections are of order one and, therefore, the 
corrections are reasonably small. 

3.2. Effective energies 

In confirming the asymptotics ©, the small- 
ness of (l/z)^-terms found in the combination 

&v / {zC mass) dcscrves particular emphasis re- 
garding the static limit computation of the b- 
mass via Fav |3I8| . since it yields an estimate of 
the error to Mb, originating from the matching 
to QCD, of only « 1% 0. The spin splitting 
LAr/Cspin (fig- vanishes for l/z ^ as ex- 
pected, exhibiting a rather large l/z-coefhcient. 
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Figure 4. Constrained fits with Cgpin to 2 loops. 

3.3. Summary 

Our successful tests of HQET show the contin- 
uum limits of the non-perturbatively renormalized 
QCD observables at finite z to meet the predic- 
tions of the effective theory. Only the functions 
C relating them to the RGIs of the latter induce 
perturbative uncertainties, but these are under 
control (except for Cgpin lacking the NNLO) and 
reveal the power corrections to dominate over per- 
turbative ones in the accessible z-range. Finally, 
our results appear promising to determine 1 / rub- 
corrections to B-physics matrix elements follow- 
ing (21 by extending the non-perturbative match- 
ing of HQET and QCD to subleading terms [51111] . 
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